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We present an Instantaneous Normal Modes analysis of liquid HF aimed to clarify the origin 
of peculiar dynamical properties which are supposed to stem from the arrangement of molecules in 
linear hydrogen-bonded network. The present study shows that this approach is a unique tool for the 
understanding of the spectral features revealed in the analysis of both single molecule and collective 
quantities. For the system under investigation we demonstrate the relevance of hydrogen-bonding 
"stretching" and fast librational motion in the interpretation of these features. 

PACS numbers: 62.60.+V, 63.20.Pw 



In recent years the analysis of Instantaneous Normal Modes (INM) of normal and supercooled liquids has given 
a sound improvement to the understanding of the microscopic processes underlying the dynamical properties of 
these systems. Applications of the method are to be found in the calculation of macroscopic quantities through 
the knowledge of the density of states (e.g. the diffusion coefficient and in the interpretation of atomic motion 
through the inspection of the eigenvectors J^-Q] . This last method appears to be a unique tool for the interpretation 
of dynamical features (single molecule or collective) of a disordered system in terms of correlated motions of its 
constituents. In particular it is interesting to explore how the presence of locally ordered units is reflected in the 
time behaviour of, for example, the velocity autocorrelation function (VACF). Attempts in this direction have to be 
found in the analyses of: i) the correlation function of the projection of the centre of mass (CoM) velocity of water 
molecules along the directions of the normal coordinates of a cluster of three molecules ii) the projection of INM 
eigenvectors onto the totally symmetric displacement coordinates of ZnCl^ - tetrahedral units Q|. Moreover the INM 
approach has been exploited to describe the dynamical features of molecular liquids (e.g. diatomic Lennard- Jones 
||, CS2 0) including hydrogen bonded systems like water |]||. On the other hand recent molecular dynamics (MD) 
simulations of HF have revealed peculiar dynamical features, e.g. a peak at ~ 50 ps _1 in the spectra of both collective 
and single molecule correlation functions ^[H^- In the collective longitudinal and transverse current spectra this 
mode appears to have an optical-like character, since its frequency is found to be independent of the wavevector. 
For the interpretation of the CoM VACF spectrum this peak has been related to the relative motion of two nearest 
neighbour molecules fl(i[| . Since nearest neighbours are also hydrogen bonded it is tempting to assign this dynamical 
feature to hydrogen bonding "stretching" . As a matter of fact HF molecules have been demonstrated to form irregular 
zig-zag chains of different size, being this peculiar clustering favoured by the geometry of the molecule and the strong 



electrostatic interaction 11 1. A dynamical characterisation of these ordered units has not yet been given, so that any 
assignment of the spectral features to particular dynamical process remains speculative. 

In the present Letter we report the results of an INM analysis of liquid HF aimed to give an answer to the above 
questions. It will be shown to which extent the presence of irregular chains is reflected in the INM spectra and how the 
appearance of optical-like modes can be understood from the analysis of the short time dynamics naturally expressed 
by the INM eigenvectors. 

The use of the INM analysis starts from the solution of the eigenvalue problem 

uj 2 Te = Ke (1) 

where T is the mass matrix of the system and K is the matrix of the second derivatives of the potential energy. Here 
e represents a multidimensional vector which specifies the instantaneous configuration of the system, i.e. in general 
the three CoM coordinates and Euler angles of the assumed rigid molecule. An average over many independent 
configurations has to be performed in order to obtain the distribution of eigcnfrcqucncies and any dynamical quantity 
derived from the knowledge of the eigenvector e. In order to avoid spurious effects originating in the fact that the 
mass matrix T, which depends on the sine of the polar angle, is not guaranteed to be strictly positive definite at each 
time, we have used the orientational coordinates described in Ref. [^2) to evaluate the derivatives of the potential 
energy. 

We have calculated the INM for liquid HF (at T = 203 K and p = 1.178 g/cm 3 ), using 60 configurations separated by 
2 ps with N = 108 molecules. The potential model, reported in Ref. Jr|, accounts for the interaction of three fractional 
charges plus a Lennard- Jones contribution and it has been shown to satisfactorily reproduce the thermodynamics and 
structure of the liquid. Due to the complexity of the interaction potential we have evaluated the matrix of the second 



1 



derivatives numerically by displacing each generalized coordinate with an increment h = 10~ 5 A since it gives a good 
stability of the matrix elements and no unusually high eigenvalues. We also checked the accuracy of the results by 
the presence of three eigenmodes with zero frequency. Our value of the increment is in accordance with the one used, 
for example, in Ref. The generalized eigenvalue problem was resolved using standard methods [jl5|. 
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FIG. 1. Spectrum of the INM, with imaginary frequencies reported in the negative part, as usual. The thick line is the total 
spectrum, and has been divided in translational and rotational components, the latter of which appears only at high frequency. 
The dotted and dashed lines correspond to the contribution parallel and perpendicular to the molecular axis respectively. 

The INM spectrum is presented in Fig. [I]. A comparison with the results for water at room temperature (reported 
in Ref. gj) points out some important features. The contribution of imaginary modes (conventionally reported on 
the negative axis) is substantially larger in HF than in H 2 (16% against 6%), a result consistent with the relatively 
higher diffusion coefficient found in HF ||Io|| . Moreover there is a much more definite separation between translational 
and rotational contributions defined by: 

Pt (uj) = /j2 E «j 2 <^-c Q )\ (2) 

\ i ^=1,2,3 / 

p«h = (E E(0( w - w ») ( 3 ) 

\ i M=4,5 / 

where ef^ is the component of the eigenvector corresponding to the frequency a, referred to the i-th molecule and to 
coordinate component fi (fi = 1,2,3 CoM, fi — 4,5 rotational coordinates). 

The rotational spectrum extends from 50 to 250 ps _1 , whereas in water it goes from to 180 ps _1 . The translational 
spectrum of HF shows a clear second maximum at ~ 50 ps _1 absent in the spectrum of water where one can only 
notice an asymmetry of the low frequency maximum with a larger content at higher frequency. We believe (and 
demonstrate in the later) that this secondary maximum is related to the "stretching" of hydrogen bond between two 
HF molecules. A comparison with the spectrum of the CoM VACF points out that the peak at ~ 50 ps _1 is 
much less pronounced and separated in the translational component of the INM. 
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FIG. 2. The ratio between p±(cj) and p\\(ui) (see Eq. ^). It should be 2 for an isotropic system. 

The strong anisotropy of the single molecule dynamics is shown by performing a projection of the translational 
component of the eigenvector along the directions parallel and perpendicular to the molecular axis: 

piiM = (£ £ K^fsiu,-^)) (4) 

\ i (11=1,2,3 / 

where is a unit vector along the symmetry axis of molecule i and we have also defined p±(u>) = pr{oS) — pi|(u>). If 
the molecular motion were isotropic, the ratio between the perpendicular and parallel contributions would be equal 
to two. Fig. H shows that below 20 ps -1 this ratio becomes larger than two going over four for modes at negative 
frequencies: a result which points out that molecules can diffuse more freely in the direction perpendicular to their 
axis rather than in the parallel one. Beyond 50 ps -1 the ratio becomes equal or less than one, thus revealing that 
the motion of the molecules in this frequency range occurs, on a large extent, in the direction of the molecular axis. 
Such an observation will be relevant in discussing the arrangement of molecule along irregular chains as revealed by 
the subsequent analysis of the INM eigenvectors. As a final remark we wish to point out that the clear separation 
between translational and rotational components can be considered as a print of hydrogen bonded systems. In fact it 
is present in water and HF but not in a linear Lennard-Jones diatomic which has been investigated in || . 

To analyse in more detail the structure of the INM we need to determine which molecules participate in each given 
mode. We assume the following criterion: the molecule i said to belong to mode a if the condition J^e^J 2 > jj 
is fulfilled. The distribution of the number of molecules per mode is reported in Fig. || (solid line) as a function of 
frequency. Knowing the particles participating in a given mode, we can determine the spatial localisation of the mode. 
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FIG. 3. Solid line (left scale): Frequency distribution of the particles participating in the modes. Dashed line (right scale): 
the distribution of the mode radii, normalized with the simulation box length (14.504 A). 

If Yi is the CoM position of particle i and S a is the set of particles involved in the mode we can define a "radius" for 
the mode as the maximum distance between two particles i.e. R a = maxjj e s Q |r^ — ry| , the distribution of which is 
reported in Fig. || (dashed line). A comparison of the behaviour of the two quantities indicates a clear correspondence 
between participation ratio and extension of the modes. A lower number of participating molecules is accompanied 
by a smaller radius of the mode. Surprisingly enough, however, in the range of ~ 50 ps" 1 the modes has an extension 
of about half of the box length even if the number of molecules participating is smaller than ten. 

A better understanding of this localisation problem can be obtained by performing a projection of a mode onto the 
hydrogen bonded chains present in the system, where hydrogen bonding is defined by the same energetic criterion 
adopted in Ref. |b|. If C c denotes the set of molecules belonging to the chain c we define the projection of the mode 
a on the chain as 



which is one if the mode a is localised on the chain c, and is zero if the mode a involves molecules not belonging 
to the chain. Clearly if we consider the maximum of this projection taken on the set of all the chains, we can see 
whether the modes are localised on some chain or not. 

The result presented in Fig. |i] (dashed line), indicates that only in some frequency range the modes show a high 
degree of localisation on a single chain, in particular in the same range where a minimum partecipation ratio occurs 
and the radius of the mode is lower (see Fig. ||) 
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FIG. 4. Solid line (left scale): Frequency distribution for the number of chains involved in the modes. Dashed line (right 
scale): Frequency dependence of the mean value of the maximum projection of the eigenmodes onto a chain. 



The idea that modes at particular frequencies are strongly correlated to the presence of chains is confirmed by 
looking at the distribution of the number of chains involved in a mode as reported in Fig. ^ (solid line) . In the range 
~ 50 ps _1 (and in the region of imaginary frequencies) we find the modes are spread over two distinct chains; at the 
highest frequencies only one chain per mode is involved. Since in this range only four molecules are participating to 
the mode (see Fig. |^) we can conclude that these modes are confined on the ring chains (tetramers) which are found 
to be particularly stable jljj. 

Having said that, it is also evident that most of the modes are far from being localized over a single chain, e.g. 
all the modes whose maximum projection is less than 0.80. This result is in accordance with other INM studies of 
network forming systems Q , where it is shown that the modes do not typically reproduce the behaviour that could 
be expected on the basis of the group properties of the network, but generally have some sort of mixed character. 

In order to characterize the spatial correlation of the molecular displacements through the INM eigenvectors, we 
have examined the following quantity: 



\ 3N ^ n(R) ^ ^ 



(6) 



where R denotes a spatial range (e.g. first or second shell of nearest neighbours, defined through the minima of the 
CoM pair correlation) and n(R) the number of pairs present in that range for a given configuration. This quantity 
is a sort of mean value of the scalar product of the displacement of the particles being nearest or next-to-nearest 
neighbours in a given normal mode and at least at short time characterises in an exact way the motion of the 
molecules participating to the particular mode. It is of course positive if the particles move in phase and negative for 
out-of-phase displacements. The results are shown in Fig. |^. 
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FIG. 5. The value of the spatial correlation between particles in various normal mode (see Eq. []). The solid line shows 
correlation between the molecules being in the first shell, and the dashed line shows correlation between molecules in the second 
shell (multiplied by a factor of 5) . The zero-frequency modes (corresponding to an overall translation of the system) have been 
dropped from the calculation. 

We notice that the nearest neighbours have opposite phases in a broad frequency range around 50 ps _1 , a behaviour 
in accordance with the presence of the "stretching" mode similar to the one observed in the solid jl6| and consistent 
with the fact that, in this range of frequency, the molecules move preferentially along the direction of the molecular 
axis. We notice that the next-to-nearest neighbours are not very much correlated in this range, a signature of the fact 
that those modes are somehow localised over a chain (see Fig. |J) and do not involve many molecules, in agreement 
with the value of the participation ratio (see Fig. ^). The present result unambiguously confirms the optical-like 
character of the mode at 50 ps _1 present in the longitudinal and transverse current spectra reported in Ref. ||. 

In a previous investigation of the collective properties of liquid HF Q it has been shown that the longitudinal 
spectra have a peak, w max , at low frequency, which changes linearly with the wavevector k. Its value is found to 
remain lower that 10 ps _1 . The corresponding phase velocity v p h — uj max /k turns out to be somewhat higher than 
the ultrasonic (hydrodynamic) counterpart, but compatible with the presence of a positive anomalous dispersion as 
in the case of monatomic liquids (e.g. liquid metals). This feature is normally interpreted in terms of overdamped 
acoustic modes propagating in the system at wavevectors well beyond the hydrodynamic range. The results of the 
present INM analysis are consistent with such an interpretation, since they show that in the range below 10 ps _1 : 1) 
the partecipation ratio is large (see Fig. ||), 2) the short time displacements of nearest and next-to- nearest neighbors 
are in phase as one would expect from molecules participating to a collective (acoustic-like) motion (see Fig. g). 

In conclusion the INM analysis of HF has revealed how the dynamical features are affected by the presence of 
topological chains of hydrogen-bonded molecules. In certain ranges of frequencies there is a strict correspondence 
between INM modes and chains, where these modes are found to be localised. This result has allowed to give a 
sound interpretation to a feature present both in the VACF and collective currents, namely a peak at 50 ps _1 in the 
corresponding spectra. Such a characteristic can in fact be assigned to the "stretching" of hydrogen bonding of first 
neighbouring molecules. This finding reveals the "optical" character of the collective mode present in the currents. 
Low frequency modes are found to be spread over the whole system involving several different chains, however they are 
confined to frequencies not much higher than 10 ps _1 in agreement with previous findings derived from the analysis of 
the longitudinal and transverse currents. The very high frequency (rotational) modes (~ 200 ps _1 ) are demonstrated 
to be confined over the highly stable tetramer chains. 

Finally, we wish to stress the fact that our present study has the potential to pave the way for an unambiguous 
interpretation of the dynamical feature of other hydrogen bonded systems (e.g. water) which are still a matter of a 
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large debate. 
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